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Abstract 

A two-dimensional (2D) assembly of nonin- 
teracting, temperature-dependent, pre-formed 
Cooper pairs in chemical/thermal equilibrium 
with unpaired fermions is examined in a binary 
boson-fermion statistical model as the Bose- 
Einstein condensation (BEC) singularity tem- 
perature Tc is approached from above. Com- 
pared with BCS theory (which is not a BEC 
theory) substantially higher T^'s are obtained 
without any adjustable parameters, that fall 
roughly within the range of empirical T^s for 
quasi-2D cuprate superconductors. 

A possible interpretation of the "pseudo- 
gap" observed in some superconductors above 
Tc is that it arises simultaneously with the for- 
mation of "pre-formed" Cooper pairs (CPs). 
We propose here that such objects emerge nat- 
urally as the nonzero-total (or, -center-of-mass) 
momentum (CMM) CPs that are entirely ne- 
glected in ordinary BCS theory. 

Consider a 2D system of N fermions of mass 
m confined in a square of area and in- 
teracting pairwise via the BCS model inter- 
action Vk^k' = —V when fJ-{T) — Hlud < 
ek-i^{= h^kf/2m) and < fi{T) + huj^, 

and zero otherwise, where k = ^(ki — k2) 
is the relative wavevector of the two parti- 
cles; Vk,k' the 2D double Fourier integral of 
the underlying non-local interaction V{r,r') in 
the relative coordinate r = ri — r2; f^{T) 
the ideal Fermi gas (IFG) chemical potential 



which at T = becomes the Fermi energy 
Ep = h'^kp/2m with kp the Fermi wavenum- 

ber; 2hojp) = h'^k'jy/m the energy width of the 
annulus centered around the Fermi circle and 
where the interaction is nonzero, with ujd the 
Debye frequency. For V > this model in- 
teraction mimics the net effect of an attractive 
electron-phonon interaction overwhelming the 
interfermion Coulomb repulsions. 

If hK = n{ki + k2) is the CMM of a CP, 
let Ek be its total energy (besides the CP rest- 
mass energy). The original eigenvalue CP |j^] 
equation for a pair of fermions at T = im- 
mersed in a background of — 2 inert, specta- 
tor fermions lying within a (sharp) Fermi cir- 
cular perimeter of radius kp, is then 



'e{ki - kp) e{k2 - kp) 

2ek + n'^K^/'im- Ek' 



(1) 



where 9{x) is the Heaviside unit step func- 
tion, and the prime on the summation sign de- 
notes the conditions ki^2 = it k| < {kp + 
k'n)^^'^ ensuring that our pair of fermions above 
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the Fermi "surface" cease interacting beyond 
the annulus of energy width 2hujp), thereby re- 
stricting the summation over k for a given fixed 
K. Without these restrictions (|l|) would just be 
the Schrodinger equation in momentum space 
for the pair. Setting Ek = 2Ep — Ak, the 
pair is bound if Ak > 0, and (|l|) becomes an 
eigenvalue equation for the (positive) pair bind- 
ing energy Ak- Our Ak and Aq should not 
be confused with the BCS energy gap A(r) at 
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r = 0. Let A = giEF)V > be the usual 
BCS dimensionless coupling constant. Here 
g{Ep) is the electronic density-of-states (for 
each spin) at the Fermi surface in the nor- 
mal (i.e., interactionless) state, which in 2D 
is g{e) = L^m/27r7i^ = g, a constant. The 
Cooper equation (|l]) for the unknown quan- 
tity Ak was analyzed in Ref. B. For zero 
CMM, i^T = 0, it becomes a single elementary 
integral, with the familiar ^ solution Aq = 

2fiuj£)/{e'^^^ — 1) valid for all coupling A. For 
small K one determines |^ for weak coupling, 
A ^ 0, that 



Ak — > Aq - -hvpK + 0{K^ 

K^o vr 



where vp = \/2Ep/m is the Fermi velocity. 
This linear dispersion relation is the 2D ana- 
log of the 3D result discussed as far back as 
1964 in Ref. Jl, P- 33 (see also, Ref. Q, 
p. 336, and |5|[) but with the 2D coefficient 
2/iT in (^) replaced by 1/2. Though commonly 
confused with the Anderson-Bogoliubov-Higgs 
(ordinary) sound mode, the linear-dispersion 
result (H) corresponds to real, moving (pre- 
formed) CPs and is distinct from the zero- 
coupling ABH (indeed, IFG) phonons de- 
scribed by hvFK/\/^. A general many-body 
formalism unambiguously exhibiting p] this 
distinction involves solving the Bethe-Salpeter 
equation for Cooper pairing based not on the 
IFG, as above, but on the BCS ground-state in 
a Green's functions scheme allowing holes on a 
par with particles. 

For Nb ordinary bosons of mass niB and en- 
ergy £k = Cs with s > and Cs a constant, 
a ternperature singularity appears at Tc 7^ for 
any Q dimension d > s in the number equa- 
tion Nb = XiKie^^'^"''^^/*''^'^ - 1]"-^ at vanish- 
ing bosonic chemical potential < when 
the number of K = bosons just ceases to be 
negligible upon cooling. It is given by 



Cs 



sT{d/2) {2TTf UB 

27r'^/2r(d/s)5rf/.(l) 



s/d 



(3) 



with ub = Nb/L'^ the boson particle density, 
and gd/s{^) the usual Bose integrals expressible 
as the series 
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9M) = T.l^ 



1=1 



C(cx), 



(4) 



where C(c) the Riemann zeta function of order 
a. The last identification in (|4|) holds when a > 



1 for which C('^) < while the series 5o-(l) 
diverges for a < 1, thus giving Tc = for d <2. 
For s = 2 and d = 3 one has C(3/2) ~ 2.612, 
and since C2 = /2mB (3) then reduces to 

the familiar formula Tc ~ S.Slh'^n^J^ /niBkB of 
"ordinary" Bose-Einstein condensation (BEC). 
On the other hand, for either particle or hole 
bosons with (positive) excitation energy ek = 
Aq — Ak given asymptotically by the linear 
term in (^) for all K, we have Ci = a{d)fiVF 
where [|] a{d) = 1, 2/7r and 1/2 for d = I, 2 
and 3, respectively. Now Tc is nonzero for all 
d > 1 — which is precisely the dimensionality 
range of all known superconductors including 
the quasi- ID organo-metallic (Bechgaard) salts 

The number of bosons in the boson-fermion 
mixture in chemical/thermal equilibrium turn 
out |l^ to be temperature-dependent, and it 
is in conserving the fermion number that the 
singularity arises. As in the case of the pure 
boson gas (3), a linear rather than a quadratic 
dispersion relation is needed to obtain BEC 
in exactly 2D. All this emerges in a statisti- 
cal model for the ideal binary gas mixture of 
bosons (the CPs) and unpaired (both pairable 
and unp airable) fermions in chemical equilib- 
rium [^2| . Thermal pair-breaking of the bosons 
into unpaired pairable fermions is explicitly al- 
lowed. At any T the total number of fermions 
in 2D is TV = L'^k], /2t: = Ni + N2, and is just 
the number of non-interacting (i.e., unpairable) 
fermions A^i plus the number of pairable ones 
The unpairable fermions obey the usual 
Fermi-Dirac (FD) distribution with the IFG 
chemical potential fi. On the other hand, the 
N2 pairable fermions are simply those in the 
interaction shell of energy width 2tujJD so that 



9{e) 



N2 = 2 de-r, ^ 



2ghujD, (5) 



where /3 = {kBT)~^, since a constant g{e) ren- 
ders the remaining integral exact. At any in- 
terfermionic coupling and temperature these 
fermions form an ideal mixture of pairable but 
unpaired fermions plus CPs that are created 
near the single- fermion energy //(T), with bind- 
ing energy Ak{T) > and total energy 



Ek{T) = 2^(T) - Ak{T). 



(6) 



This generalizes the T = definition Ek = 
2Ef — Ak given below (1). 

The Helmholtz free energy F = E — 

TS, where E is the internal energy and S 
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the entropy, for this binary gas "compos- 
ite boson/pairable-but-unpaired-fermion sys- 
tem" at T <Tc is then readily constructed [11| 
in terms of: a) n2(e), the average number of 
unpaired but pairable fermions with energy e; 
b) NsfliT), the number of (bosonic) CPs with 
zero CMM at temperature T; and c) Nb,k{T), 
that number of excited pre- formed CPs (i.e., 
with arbitrary nonzero CMM K) and a cutoff 
Kq physically defined Q by Akq = denot- 
ing the value of K beyond which a CP breaks 
up. The free energy F2 of the pairable fermions 
is to be minimized subject to the constraint 
that ^^^^2 is conserved. If N2o{T) is the number 
of pairable but unpaired fermions, the relevant 
number equation for the pairable (i.e., active) 
fermions is then 

N2 = N2o{T) + 2[NB,o{T) + NBfi<K<Ko(T)] 
= N2o{T) + 2Nb{T), (7) 

where NBfi<K<Ko (T) denotes the total number 
of "excited" bosonic pairs (namely with CMM 
such that < K < Kq), i.e., NBfl<K<Ko{T) = 
Yl,o<K<Ko ^b,k(T). At r = two distinct cou- 
pling regimes emerge: a) for Aq < 2huJB) or 
for A < 2/ln2 ~ 2.89, we have that iV2o(0) = 
g{2TLUj£) — Aq); while b) for Aq > 2Tlujd (or 
A > 2.89), iV2o(0) is identically zero. Hence, 
the number of bosons Nb{0) at T = from 
(7) is just Nb{0) = l[N2- A^2o(0)]. Using 

ffor N2 the fractional number of pairable 
mions that are actually paired at T = 0, 
namely 2Nb{0)/N2 = 1 - A^2o(0)/iV2, becomes 

simply Ao/2nuJD = {e^^^ - 1)"^ — > 6"^/^ as 
A ^ 0, for A < 2/ In 2 ~ 2.89, and unity for 

A > 2/ In 2. As Nb{0) = ^Aq for A < 2.89, 
only those fermions in an energy shell of width 
^Aq around the Fermi surface actually pair at 
T = 0, while for A > 2.89 all pairable fermions 
pair up since then Nb{0) = gfiuJo = ^^2- For 
T > 0, 2Nb{T)/N2 = 1 - N2o{T)/N2 de- 
creases with T, provided one knows Aq{T) for 
any T > and ascertains that it decreases. For 
r > 0, the e{ki - kp) = e{ek^ - Ep) in (|) be- 
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comes 1— n(^fcj, where n(^fcj = (6^^*= 
the FD distribution with = e^^ — /u(T), with 
the IFG chemical potential /i(r) in 2D given 
exactly by ^{T) = (3'^ ln(e^^^ - 1) — > Ep as 
T — > 0. Similarly for 0{k2 — kp). Since ki = k2 
implies that = ) (S) then leads to a sim- 
ple generalization to nonzero T of the K = 
CP equation, 

1 = A / de(e-^« + l)-2[2e + Ao(r)]-i. (8) 



Numerical solution shows Aq(T) to indeed be 
monotonic-decreasing in T for any fixed A and 
hujD- Further, the solution of Aq{T*) = is, by 
inspection, T* = 00; this infinite "de-pairing" 
temperature is unrealistic and undoubtedly an 
artifact of the simplest version (1) of Cooper 
pairing used here as a starting point; see, how- 
ever, Ref. 1^. 

Modeling our system as a pure boson gas of 
CPs (i.e., neglecting the background unpaired 
fermions) but with a temperature-dependent 
number density n b (T) converts the explicit Tc- 
formula (3) into an implicit one. For s = 1 and 
d = 2 it becomes, since (72(1) = C(2) = tt^/G, 



Tr 



4^/3 hvp 
7r3/2 kB 



UBiTc). 



(9) 



This requires UBiT) = Nb(T)/L'^ which in 
turn requires Aq(T) as determined from (^), 
and follows from the expression 2Nb{T)/N2 = 
1 — N2o{T)/N2. Solving this self-consistently 
with (7) for A = 1/2 gives the remarkably 
constant value Tc/Tp ~ 0.004, where Tp = 
Ep/kB, over the entire range of = huJo/Ep 
values 0.03 — 0.07 typical of cuprate supercon- 
ductors. On the other hand, the BCS for- 
mula Tf^^ ~ 1.13GDe-^/^ with A = 1/2 gives 
Tc/Tp = 0.005 to 0.011 over the same range of v 
values. Obviously, both sets of predictions are 
too small compared with the empirical cuprate 
range Tc/Tp ~ 0.03 - 0.09 [|l3|. 

The exact Tc without neglecting the back- 
ground unpaired fermions requires the ex- 
act CP excitation energy dispersion relation 
£k{T) = Aq(T) — Ak(T) which is neither pre- 
cisely linear in K nor independent of T. To 
determine AxiT) we need a working equation 
that generalizes Ref. Q for T > via the new 
CP eigenvalue equation (|8|). At T = Tc both 
Nb,o{Tc) — and ^b{Tc) — so that one gets 
[11 1 the implicit Tc-equation for the binary gas 
mixture 



V 



1 + e 



1 + g-{Ao(rc)/2+i/}/rc 
dn- 



+ 



g[Ao(T,)~A4T,)]/Tc _ 1' 



8(1 + z^) 



(10) 



where quantities with tildes are in units Ep 

or Tp] K = K/2{kj, + k^)^ with ko defined 

through huJD = h?k'jy/2m; and v = Qo/Tp. To 
obtain Tc from the finite-T dispersion relation 
one must numerically solve four equations self- 
consistently for each A and v, namely (10) in 
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conjunction with (8) for Aq{T), and Eq. (35) 

of Ref. for both Ak(T) and the breakup 

value kq{Tc). For A = 1/2 and the range of 
values 0.03 — 0.07 typical of cuprates, the re- 
sulting Tc/Tp falls within the aforementioned 
empirical range 0.03 — 0.09 [|l3l|. For cuprates 
d ~ 2.03 has been suggested |I14|| as more real- 
istic since it reflects inter-CuCMayer couplings, 
but our results in that case would be very close 
to those for d = 2 since, e.g., from (3) Tc for 
s = 1 (but not for 2) varies little with d around 
d = 2. In fact, if mB± and niB are the boson 
masses perpendicular and parallel, respectively, 
to the cuprate planes, an "anisotropy ratio" 
t^B/TUBi. varied from to 1 allows "tuning" 
d continuously from 2 to 3. 

Other boson-fermion models |15|[16|[17|[18| 
I p^ have been introduced, some even address- 
ing |17]||l8[| d-wave interaction effects as op- 
poseoto the pure s-wave considered here, and 
some also focus [18|[|^] on the pseudogap. But 
calculating cuprate Tc values in quasi-2D with- 
out adjustable parameters has not been at- 
tempted, and indeed Tc = is predicted in 
exactly 2D. 

In summary, a simple statistical model 
treating pre-formed CPs with both zero and 
nonzero CMM as non-interacting bosons in 
chemical/thermal equilibrium with unpaired 
fermions is proposed that gives rise to a bo- 
son number that is strongly coupling- and 
temperature-dependent. Since the CP dis- 
persion relation is approximately linear for 
nonzero CMM, it exhibits a BEC of zero-CMM 
pairs at precisely 2D. In contrast to both BCS 
theory — which is not |2C] a BEC theory — and 
simpler BE models excluding either the break- 
able character of the pre-formed CPs or the 
presence of unpaired fermions, exact Tc's for 
the boson-fermion binary mixture based upon 
the exact CP dispersion relation are found in 
rough agreement with empirical cuprate values 
with no adjustable parameters. 
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